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1 Introduction 


Higher loop calculations present great difficulties especially in the presence of masses. When 
there is a large number of parameters (masses and external momenta) in the problem, analytic 
expressions can hardly be obtained beyond one-loop. In such cases one has to resort to numerical 
or approximative methods. In particular the rather large spectrum of masses in the Standard 
Model requires numerous and different approximation schemes (e.g. small fermion masses 
mj/m| -C 1, close vector boson masses (m| — -C 1, large top mass m|/mj -C 1 etc.) 

and often a problem is reduced to diagrams with one mass. 

In this work we consider the evaluation of 2-loop diagrams of self-energy and vertex type 
with one non-zero mass and all others zero. 

For the analytic evaluation of ‘master integrals’ with masses various techniques have been 
developed. Such calculations in QED are known since [|l|. Using dispersion relations 2-loop 
2-point integrals relevant for QED and QCD were obtained in |^. Results for a class of 2-loop 
self energies are given in 0-0. The method of differential equations was developed in 0]. 
Some recent results, applying this method for 2-loop vertex functions, are presented in [^. At 
the 3-loop order analytic results are known only for vacuum bubbles so far [^. 

Very powerful approaches for the evaluation of diagrams, which are in fact our starting point, 
are the asymptotic expansion |]^ and/or Taylor expansion [^, in an external momentum 
squared, which are identical if only one non-zero mass is involved. These methods allow to 
expand a diagram in a series with coefficients, which in many cases can be represented in a 
relatively simple analytic form. Analytic continuation by means of a mapping in the plane of 
momentum squared allows to calculate the diagrams on their cut. For propagator type diagrams 
such expansion was successfully applied in |^, |T^ while for vertex functions it is discussed in 
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So far the method of expansion is considered as semianalytic in the sence that only a 
limited number of coefficients can be obtained explicitely. In this paper, however, we want to 
go one step further. We calculate the hrst few coefficients of the expansion of a diagram with 
only one non-zero mass by means of the asymptotic expansion m method. The method of 
differential equations 0 then yields an idea, like in |^, what the general analytic form of these 
coefficients might be, providing some ‘basis’ in terms of which they might be expressed. The 
Ansatz of equating the explicit coefficients obtained from the asymptotic expansion to a linear 
combination in terms of the basis elements, yields a system of linear equations, which can be 
solved to yield the desired representation of the coefficients. 

The main problem in this approach is the choice of the basis elements. We start from 
so-called harmonic sums which are particularly relevant for moments of structure functions in 
QCD (e.g. [^). These functions are directly related to (generalized) polylogarithms [|T^ and 
therefore it is not surprising that they appear in the analysis of massive diagrams (see also 
recent research about these sums 0)- However, not every massive diagram even of self-energy 
type can be expressed in terms of these sums. In particular, if a diagram possesses a cut with 
three massive particles in the intermediate state then it is rather a subject of elliptic integrals 
than polylogarithmic ones P, |T^. The method should work in this case as well provided one 
hnds the relevant class of basis elements. However, in this paper we do not consider diagrams 
of such type and concentrate ourselves on those with one- and/or two-massive-particle (m- and 
2m-) cuts. We give the results for all selfenergies and a certain number of vertices with m- and 
2m-cuts. In particular we include diagrams relevant for Z-decay. 


1 
















The paper is organized as follows: after introducing some notation in Sect. 2, we explain 
in Sect. 3 the main idea of our approach in terms of some simple examples. Sect. 4 contains 
a ‘straightforward’ prescription for the generalization of our basis elements, while in Sect.5 
we demonstrate how to hnd more complicated ones occurring in 2m- cuts. The differential 
equation method (DEM) in Sect. 6 demonstrates a further method to explore the form of basis 
elements. Sect. 7 contains a proposal of how to perform the large g^-expansion by means of the 
Sommerfeld-Watson transformation, providing an explicit example. The Appendices contain 
further notations and results. 


2 Notation and definitions 


The generic topologies considered in this work are displayed in Fig. 1. Lines of the diagrams 
represent scalar propagators {p^ — m? + with either mass m or zero mass. For both 

ultraviolet and infrared divergences we use dimensional regularization in the Minkowski space 
with dimension D = A — 2e. Each loop integration is normalized as follows 

/ dk = {m\^^f j ^ (1) 


with 7 ^ the Euler constant. 

Our starting point is the large mass expansion |]^ of the diagrams. Here we just note that 
the result of a large mass expansion for these diagrams reads [z = /m?) 


J 




j=0 


k=0 


( 2 ) 


where a is the dimensionality of the diagram, oj and u independent of n are the highest degree 
of divergence and the highest power of log(—respectively (in our cases a;, i/ < 4). The 
coefficients j ^ are of the form ri + ( 2 'f '2 + ... + Cu'f'i, with being rational numbers and 
Cc = C(c) is the Riemann (^-function. 

Series (^ always has a nonzero radius of convergence, which is defined by the position of 
the nearest nonzero threshold in the g^-channel. For brevity we shall call m-cuts (2m-cuts, 
etc.) possible cuts of a diagram in g^ corresponding to 1 (2, etc.) massive particles in the 
intermediate state. 


3 Diagrams with m-cuts 


We start from diagrams with the simplest threshold structure, i.e. with m-cuts and (possible) 
0-cuts. Such diagrams are particulary easy to handle and results for all of them within the 
considered generic topologies of Fig. 1 can be found in the appendices. It will be shown that 
all are expressible in terms of harmonic sums. 

As a first example consider the 2-loop 2-point function of Fig. 2, 7 i 3, which was considered 
in [^]. Using the standard large mass expansion technique one can get the first few coefficients 
of the expansion of this diagram in powers of ^ = g^/m^ 


— anZ'"- 
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where (2 = C(2) is the Riemann (^-function. 

It takes seconds to evaluate the first 10 terms (we use a package written in FORM p3| ) 
but the evaluation time increases exponentially for the higher order terms. Although 10-20 
coefficients are sufficient to get a very precise numerical value for the diagram, even at \z\ ~ 10 ^ 
(via conformal mapping and accelerated convergence [Mj), this cannot be called a satisfying 
result. 

Our key observation is that the information contained in the first 10 (or even less!) coeffi¬ 
cients of the series (^ is enough to restore the whole series. Namely, each coefficient a„ can be 
expressed as some combination of harmonic sums and their extensions. In this case the series 
can be continued numericaly or analyticaly to any point of the complex z-plane. 

Thus our first step is to find an expression for the higher order terms of the series (^ as a 
functions of n i.e. = a{n). 

To achieve this let us search for a„ as a linear combination of elements of a certain kind. 
For these we take 1 /n and Skiji — 1) = 1/j^ (harmonic sum). As a matter of observation, 

the rule is that one has to take into account all possible products of the type CaSb-^ ■ ■ ■ Sb^^/n^ 
with the ’transcendentality level’ (TL) a -|- 61 6 ^ -|- c = 3. It is obvious that one can 

exclude ^3 beforehand since it never appears on the r.h.s. of (^. Thus we have the following 
Ansatz for 


C 2 


ttn — —Xi 4- C 2 S 1 X 2 -l- S^x^ -|- 8281 X 4 ^ 4- 8 iX^ 4- 
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>^1 >^1 1 

Xq 4- Xj 4—4—rXg , (4) 


where xi,..., xg are rational numbers independent on n and we omitted the arguments in the 
S'-functions. These can be taken with argument n or n — 1. We choose the latter option. 
However, either choice gives a solution since the difference is only in rearranging 1/n terms. 
Later on we refer to the structures in (^ as ‘basis elements’ or just ‘basis’. Indeed, the functions 
Ca 8 bi{n — 1)... S'bj,(n — l)/n'’’ are algebraically independent. 

Inserting the expression for an from (^ into the l.h.s. of (||) and equating equal powers of 
z, we obtain a system of linear equations for the Xj. One needs at least 9 first coefficients on 
the r.h.s. of (^ to solve the equations for 9 variables and any system of more than 9 equations 
should be consistent and have the same solution if the Ansatz (^ is correct. An explicit 
computation ensures that the system can be solved in terms of rational numbers for the x* and 
this latter consistency property holds. The solution is xi = 2,X6 = 2,X8 = —2, X 2 , 3 , 4 , 5 , 7,8 = 0 
i.e. the answer for the diagram at hand is 



.C 2 , . ^2(n-l) , gi(n-l) - 

n n n? 


(5) 


This is the desired result. One notices that in (^ and Ansatz (^) the terms proportional to 
(’2 are linearly independent and can be treated separately, i.e one could solve a system of linear 
equations for xi, X2 and another one for X3,..., xg. In this sence only the 7 first coefficients of 
the large mass expansion are required. 

Series (^) converges for \z\ <1 and represents in the domain of convergence an analytic 
function. It can be analytically continued into the whole z-plane except for the cut on the real 
axis starting at z = 1 where /13 has a branch point (threshold). In particular for the given 
diagram, summing the series, we find the representation 


As = ^ - 2 C 2 log(l - z) -2 log(l - z) \^i 2 {z) - OS'!, 2 ( 2 ;) 


( 6 ) 
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where is a polylogarithm and Sp^q{z) are the generalized Nielsen polylogarithms (see 

also Appendix C). 

Both representations (|^) and (@) are eqnivalent in the sense that they both represent the 
same analytic fnnction and can be nnambiguously converted into each other. However, we find 
it more convenient to use the series representation and refer for the continuation to App. C and 
E. Moreover a series representation is usually shorter than a representation in terms of known 
functions and/or an integral representation. 

Equation can be considered as expansion of the general coefficient in terms of a 
‘basis’ with (unknown, rational) coefficients x*. Now the question arises: what are the basis 
elements in general? Definitely there must be a certain connection between the topology and 
the structure of thresholds of a diagram on the one hand and the structure of the basis elements 
on the other hand. So far we are lacking rules for predicting a basis of a given diagram. The 
power of the method, however, is that given a set of basis elements for one diagram, it can 
be used to find the solution for other diagrams. Often, though, one has to ‘generalize’ already 
known (harmonic) sums in a basis. In Sect. 4 we explain how to do this. 

In particular it will be shown that many 2-point and 3-point functions have similar bases 
and thus solving some 2-point integrals (which are simpler) we can find solutions for 3-point 
integrals. As an illustration consider the diagram shown in Fig. 3, P 5 . Its large mass expansion 
looks like 


P. = 




£ 2 :” (r® log^(- 2 ;) + log(- 2 ;) + + r^n ‘^\2 + r 


( 0 . 0 ) 
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(7) 


with the r’s being rational numbers. It is obvious that one can search for a solution for each of 
the r’s independently. Again we use the same set of functions (l/n“ and Sb) as above but now 
with different transcendentality level(s). The system of equations has a solution only if we add 
the factor (—)” which can be seen easily by inspection of the series. At the end we arrive at 


P.= 


1 c» / Q c C2 C 

^ + ^ - 2%) log(-j) 

[q ) \n ^ n r). 


-6*+ 2^+ 6^ 

n n n 


c C c C ' 

n n'^, 


( 8 ) 


where as above we take all harmonic sums (S's) with the argument n — 1, which is omitted. 


The similarity of (|^) and (H) is spectacular. There are two comments in order. 

At first we point out that in (H) the part without log(—z) has basis elements (^aSb/'rf obeying 
a -|- 6 -|- c = 4 (i.e. it has a basis with TL= 4). Terms proportional to log(—z) are of 3rd level 
while those proportional to log^(—z) of the 2nd. One can say that log(— 2 ;) itself is of 1st level 
and each log“(—z) reduces the level of basis elements by a units. This is the general behaviour 
for all diagrams we consider. The same rule applies to UV and IR poles 1/e: if they are present 
in a diagram i.e. each l/e°' reduces the level of basis elements by a units. 

Secondly we wish to underline the presence of the factor (—)”■. Strictly speaking the basis 
now is not CaSb/n^ any more but rather {—)^(aSb/n'^- In principle one can expect that a mixture 
of both bases occurs and this is indeed the case for some diagrams. Moreover, the (—) may 
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stand not only in front of Sa but also inside the harmonic summation. Thus we are led to the 
alternating harmonic sum^ Ka{n — 1) = 

It will turn out that the results obtained above indicate already a general property, i.e. in 
2-loops all the selfenergy functions have TL=3 and all vertices have TL=4 (at least for the 
diagrams under consideration). 


4 Generalization of basis elements 


It was mentioned above that we cannot directly predict the basis of a diagram. Therefore we 
experiment with basis elements of diagrams calculated before to match the expansion coeffi¬ 
cients of other ones. Another extremely useful trick is that we hrst establish the structure of 
the lowest level in a diagram (i.e. the structure of the coefficients of the highest pole in e or 
the highest power of log(—; 2 )) and then try to reconstruct terms of higher level. If, say, the 
l/e’^ coefficients are expressed symbolically in terms of basis elements f^a) with TL=a then the 
part is expressed in terms of some /(a+i) etc. Therefore we need rules of generalization, 
i.e. to find higher level elements from lower level ones. These are discussed in this Section. 

We start from the following observation. Let J be a D-dimensional Feynman integral 
(properly normalized to be dimensionless) depending on z = /m?. It can be formally written 

as multiple series of the hypergeometric type 


J{z) 


r(ai{j,D})...r(a^{j,D}) 

n.w.r(/3i{j,/^})...r(/3,{j,D}) 


(9) 


where T is the Euler T-function. Here the symbols ai{j, D}, -0} yjj, D} stand for 

some linear combinations of the summation indices jk and the space-time dimension D = 4: — 2e. 
The form (^), though rather obvious by itself, can be deduced e.g. from the Feynman parameter 
representation or a-representation. From this point of view, ’to evaluate’ diagram J means to 
bring multiple sums in (^ to some ’known’ hypergeometric form of series or at least to remove 
as many summations as possible. 

Usually, however, one is interested in the expansion of a diagram around £ = 0. On expan¬ 
sion, (H) develops some singularities in the z-plane and possible poles 1/e:. The expansion of 
F-functions produces -^-functions or equivalently S'-sums. Namely for j integer 

p!'! ^ = (i - l)!exp|£a''£^5fe(j - 1)| , (10) 

r{l + ae) J 


with Sk{n) = X]s=i 1/s^. This shows that harmonic sums naturally play a special role in our 
subject. 

Terms of lowest order are obtained by replacing the exponent in ([lOD by unity and (P) 
becomes a multiple sum of factorials. Then the coefficient of z"' is a hnite sum over factorials. 
Terms of higher level are obtained by expanding the exponent in (IT) in e, which produces 
S'-functions in (j^). 

This suggests the following recipe of generalization: if /(„) is found to be a multiple (finite) 
sum f(a){ji) = jyij}C{j}, then higher level objects are obtained either by multiplying /( q ) by 

^These sums were used in 
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S'-functions or by inserting S'-fnnctions nnder the summation sign of /(a) to give 


f{a+i){n) = 

{il 

n n 

/(a+2)(n) = /(a+l+l)(n) = 

{j} {j} 

etc., 


( 11 ) 


where the arguments of Sk are some linear combinations of summation indices. The same 
applies of course to the function Sa itself giving Sa,b etc. 

From (prl) it is seen that any /(„) produces in general many functions f(a+i) and higher ones. 
In the case of a multiple summation the number of possible combinations of may get 

too large and this general rule is getting rather formal than practical. However, in practice 
we always met the situation that fa is a simple (one-fold) sum. Then, inserting Sk{j — 1) and 
Sk{j) or equivalently Sk{j — 1) and 1/j, (B becomes 

f{a+i){n) = Cj 5i(j - 1)|, 
j=i {3 ) 

/(a+ 2 )(n) = E c, 14’ Slij - 1), S,{j - 1)| , 

j=i k 3 3 j 

etc. (12) 


The choice ([l^) is justified a posteriory. 

This prescription, however, must be modified if a diagram possesses 2m-cuts (with a thresh¬ 
old at = 4m^). Such diagrams produce in the sum (P) functions like Viflj + ae) and according 
to (pT|) S'fc(2j — 1). In such cases, along with Sk{i — 1), we insert also S'fc(2j — 1). 

In conclusion of this Section we want to stress that the recipes formulated above do not 
necessarily work in every case for any diagram. But in many cases they do. A possible reason 
of a failure of the prescription (Q) may be e.g. the following : some new basis element begins 
to contribute only, say, at the order 1/e while at the order l/e"^ it happens to drop out. In 
such cases it will be missed in the analysis and one has to use other more powerful and direct 
methods of determination of low level terms. 


5 Diagrams having both m- and 2m-cuts and Vl^-functions 

We now establish the structure of the lowest level for the set of diagrams having both m- 
and 2m-cuts. While for diagrams with only m-cuts the complete basis consists of elements 
1 /n", Sb, Kc and their generalization due to ([T^) , in the presence of 2m-cuts the basis changes 
drastically. In this case the behaviour of the coefficients is governed by a new class of elements 
which we call W. The generic structure of the lowest level is 

( 13 ) 

\3 ) 3 

and higher level functions are obtained following the rules of the previous Section. 
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The hh-structure can be deduced from different diagrams with 2m-cuts but most easily we 
hnd it from the nonplanar graph shown in Fig. 3, A^i 2 . Namely this diagram has a double 
collinear pole and we can easily perform both integrations in calculating the 1/e^ contribution. 
Consider at first the massless lines 3 and 5 (see (J^) for the normalization) 


iVl2 = 


dki 

ki{ki -piY 


(j){ki,k2,PuP2), 


(14) 


where 0 stands for the rest of the diagram. Introducing the Feynman parameter ai this can be 
written as 


Nu = 



dki 

(fcf - 2aifcipi)2 



dki 

{ki - aipi)‘^ 


(15) 


where we have used pf = 0. This integral diverges in 4-dimensions when {ki — aipi) —> 0, 
i.e. there is a collinear divergence. For the calculation of the most singular term in e: one can 
effectively replace 


1 

{ki - aipi)'^ 


i—5^^\ki - aiPi), 
e 


(16) 


and the same replacement applies to the pair of lines 4 and 6 and integration dk 2 - Therefore 


iVl2 


l/e^ 



Integrating out ki, k 2 with the help 
we arrive at {z = /m^) 


dki dk2 d^^^ki — aiPi) 6^^\k2 — o>2P2) 

[{ki - k 2 - pi^ - m^] [(/c 2 - ki- P2Y - m^] ‘ 

of the (5-functions and using p\ = P2 = P1P2 


(17) 

gV2, 


Ilf 1 f dai da2 

Jo (1 — 0(1 (1 — 0 ( 2 ) 2 :)(1 — (1 — ai)Q; 22 :) 

^ . jnir .^M2j)!l 

(2^-1)! 

^ KQ) ^i J ^ 


(18) 


If we assign to the factor the 0th level, then the expression on the r.h.s of (p!8D has 

the correct 2nd level as it should be for the double pole part of a vertex function. Using the 
results of Sect. 4 we obtain functions of higher levels hF 2 , hFi^, W 3 etc. (see Appendix A). 


6 Differential equation method 

It has been demonstrated above that the most important point in our approach to set up a 
Taylor series for the diagrams under consideration is to find the proper basis elements. In the 
last Section we have seen a special calculation to find one more basis element. There is no 
general way of how to find further such elements but another effective approach is that of the 
differential equation method (DEM) |^. It was used in to find the expansion of diagrams 
Fi 26 and in the notation of Fig.3. Therefore we give in this Section a short review of 


7 














the DEM and demonstrate the Ending of another basis element by calcnlating the two-point 
diagram J 15 of Fig. 2. 

The DEM allows one to calcnlate massive diagrams by redncing them to others having an 
essentially simpler structnre More precisely, one writes a linear difierential eqnation (w.r.t. 
masses) for the diagram nnder consideration. The inhomogeneons term of the eqnation is a snm 
of diagrams with smaller nnmber of lines. These in tnrn can be rednced via another difierential 
eqnation to even simpler diagrams and so on. Solving the so obtained eqnations iteratively one 
gets for the initial diagram an s-fold integral representation which is in general simpler than 
the usnal parametric or dispersion integrals. Often some of the integrations can be performed 
analytically. 

For completeness we give some formnlae that will be nsefnl in the fnther discnssion. We 
introdnce the following graphical representation for the propagators 

_._1_ = a _._1_ — a O') 

WTW ’ (g" - + *0)“ ^ ’ 

a and m are called the index and mass of this line. Fnrther (except when mentioned otherwise) 
all solid lines have the same mass m. Lines with index 1 and mass m are not marked. 

The basic tool is the integration by part relation ||^ which is obtained by mnltiplying the 
integrand of a diagram with {d/dk^)kfj, and using / (i^div(...) = 0 . As a result for a triangle 0 
with arbitrary masses one has (here lines with index at have mass rrii) 



We stress the fact that the basic line of the triangle plays a special role (in the following we call 
it ’distinctive’). This relation is used in the following manner. All triangles on the r.h.s. with 
’shifted’ indicies ai + 1 can be written as derivatives of the initial triangle w.r.t. rrif. Generally 
speaking, one has to consider all masses ruj to be different and put them equal at the end of the 
calculation. In practice one can often eliminate some of the derivatives by combining several 
difierential equations. 

Whenever a self-energy subdiagram occurs the following relation is used 


Qfi, mi 



01.2 , ^2 


r(ai-|-tt 2 ~ D/2) D ds 

' r(ai)r(a2) io (1 - s)ai+l-D/2ga2+l-W2 


g 


a\+a2—D/'2 



( 21 ) 


This allows one to reduce an /-loop to an (/ — l)-loop diagram with one propagator having 
’mass’ m^/(l — s) -|- m\/s. 

We now apply this technique to the evaluation of diagrams with one mass. Consider first 
the symmetric planar diagram P 126 (Fig. 3). It possesses only 2m-cuts and zero-cuts but no 
m-cuts. Therefore the nearest singularity in the complex z = /m? plane is placed at z = 4 


®A similar technique has been introduced recently by Remiddi [ pl| . 

^Similar relations may easily be obtained for arbitrary n-point diagrams (see e.g. and discussion in p^). 















and the diagram obviously has some new structures different from those already discussed. 
Applying (20) with distinctive central rung we have 



It is easy to see that the last two terms on the r.h.s. can be combined into the total derivative 
of the initial diagram w.r.t. Therefore we immediately obtain the differential equation for 
our diagram. The inhomogeneous term now is the sum of the first and the second diagram on 
the r.h.s. of (22). The former is trivial while for the latter we repeat the procedure with the 
left vertical line being distinctive: 



(23) 


The resulting diagrams on the r.h.s. all have self-energy insertions which can be replaced 
according to (21) and the remaining triangles are solved by usual Feynman parameters. 

After some rather long but trivial transformations we obtain 


Pl2G — 
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+ log[s(l - s)] log2(l - 0 - ^log'(l - 0 - 12^1,2(0 - 51og(l - 0Li2(0 
+ ^ log(-^) log^(l - 0 + ^ ^og^{-z) log(l - 0 


(24) 


were ^ = zs(l — s) and *S'i^ 2 (' 2 ) is the Nielsen polylogarithm (see Appendix C). The above 
remaining integral is very difficult to evaluate analytically. However the expansion in z is 
welcome 


7^126 — 


(7 


2'|2 


E 


n=l 

3^ 15 

2 


'2rv 


n 




+ 


1 r 


£ L 


-Si - log(- 2 ;) 


-5S2 - ySf + iSiSi + 2^ - S, log(-2) + ilog'^(-..)|, 


(25) 


with Sa = 5a(n — 1 ) and Si = S'i(2n —1). This expression has the correct threshold at = 4m^ 
due to the factor . The appearence of S is somehow expected as was mentioned at the end 

of Sect. 5, but the factor is difficult to guess without additional information (compare 

however with (p! 8 D ). 

We now turn to the 2-point asymmetric diagram /15 (Fig. 2) which contains another new 
basis element. Using (20) with distinctive central rung for the left triangle we obtain the 
following differential equation 


2 
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2''-LZ 2''- 


One gains little from this because of the last term on the r.h.s. This has the structure of the 
initial diagram but with one massless propagator squared. To get rid of this term we can write 
another equation, using the second massive line as distinctive of the triangle: 


D-A- 2m^-^ 
dm 



+ 



+ 





^ 



Taking a proper linear combination of these equations, we can exclude the bad last term. In the 
resulting equation the inhomogeneous term consists only of diagrams with self-energy insertions 
which can be solved immediately. As a result for J 15 we write its expansion in z 


hb 


E 

n=l 



n 



n—1 


log(-z) - C 2 + 3 ^ 



1 2 ^ 
~ H—§■ + 4 



(26) 


The inner sum in (^) is the new type of basis elements which we discovered from this 
calculation. In a somewhat more general way of writing it is 


Vain 


n—1 

i) = E 

i=i 


^2j\ E 

. 3 ) 


(27) 


Further basis elements were obtained in a similar fashion. They are listed in the appendix. 
We do not give all the details of their derivation, which was done in a way analogous to the 
above. 


7 Analytical continuation 

Taking care of the zero-thresholds by factorizing log(—z), [z = q^/M"^) as described above, all 
series under consideration have finite radius of convergence \z\ < |^o| where zq is the lowest 
non-zero (pseudo)-threshold. Therefore these series need to be analytically continued on their 
cut since this is in general the region of physical interest. The most straigthforward way is the 
continuation via conformal mapping and special summation techniques (usually Fade approxi- 
mants) []^. In principle this allows one to continue the diagrams to the whole complex ^-plane 
including the cut; at least this method works perfectly for low and moderate \z/zq\. For the 
diagrams we consider here, the situation is even better than in cases with several masses since 
our analytic results for the Taylor coefficients allow to calculate arbitrary many of them easily 
with any desired precision. 


Nevertheless it is of interest to investigate the domain \z/zo\ S> 1 and have series in powers 
of 1/z. The first few coeficients of such series can of course be obtained by expanding any 
diagram in the ‘standard’ manner |M, |T^, [^] in large q^. This expansion is known to be 
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much more difficult than the corresponding g^/M^-expansion. One of the reasons is that the 
number of contributing subgraphs in this case is usually larger. In the language of the present 
paper the ‘basis’ is more extended (see examples below). 


Here we present a different approach for the diagrams with one non-zero mass, starting from 
the series in z = The idea is to apply the Sommerfeld-Watson transformation to the 

z-series : given J = CnZ^, it can be written as a Sommerfeld-Watson contour integral 


J 


1 r CnZ^ 
2ni J sinfTrn) 

C 1 +C 2 


(28) 


The integral in (^) is to be performed along a closed contour in the n-plane consisting of two 
pieces: Ci is a line parallel to the imaginary axis (7 — ioo, 7 -|- ioo) with 0 < 7 < 1 and C 2 is a 
semicircle closed in the right-half-plane at inhnity. 

The representation (^ 8 ]) is valid of course only if the coefficients c„ (considered as analytic 
function of n) decrease fast enough at inhnity and do not have singularities to the right of Ci. 
Moving C\ to the left in the negative direction, one has to add the corresponding residues. 
We now suppose that c„ is regular in the whole n-plane except (possible) poles at nonpositive 
integers n = —A:, fc = 0,1, 2,... and has the required behaviour at inhnity. Then 


J 


00 


E 


res„=_fc 


Aixu ^ yU 

sin(7rn) 


(29) 


The crucial point now is to show how to continue analytically each of the structures occuring 
in this paper. First let us note without proof that all singularities of these functions are placed 
on the real axis at nonpositive integers. Therefore we need to continue them through the 
neigborhood of these points i.e n = —k -|- 5 with /c = 0 , 1 , 2 ,... and |5| < 1. 

To start with, let us consider one-fold sums. First of all we note that Si is continued by 


5'i(n - 1) = V’(n)-f-7^; (30) 

where 'ip{n) is the logarithmic derivative of the F-function. All sums Sa with a > 1 are continued 
via the ■^-function and its derivatives. All Ka are simply expressed through Sa 

Ka{n) = Sa{n)-2^--Sa{[n/2]). (31) 

To continue the other functions we use the following general method. Let /(n — 1) = 0(j)- 

Then, formally 

C30 00 00 

fin - 1 ) = E - E ^ij) = fi^) - E + n)- (32) 

j=l j=n j=0 

By inspection of the elements under consideration we see that 4>{n) is always an analytic 
function in the whole n-plane apart from (possible) poles at nonpositive integers. Then, with 
n = —k -|- S, we obtain for the last term in (^) 


00 A:—1 00 

E 0(j - + ^) = E ^if -k + s) + (j){5) + Y. 

j=0 j=0 j=k+l 

k 00 

= E f^i-f + <^) + (fi^) + E f^if + '^)- 

j=i j=i 


(33) 
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It is convenient to introduce the notation — 1) = for fhe sum with ’shifted’ 

summation index. Combining (^) and (|33D we hnaly get 

k 

f{-k-l + 5) = f{oo) - f^^\oo) + 5) - ct){5), (34) 

i=i 

with k = 0,1,2,.... The remaining hnite sum can be easily transformed in our cases to a 
’known’ function and (possible) poles in 6 can be explicitely separated. In this way we obtain 


the following folmulae 

-! + «) = S„(oo) - SP(oo) - (-rs>-‘>(k) - i (35) 

K4-k - 1 + i) = ii'.(oo) - (36) 

V^(-k -! + «) = 14(«d) - VW(oo) + (-r^ctg{wS) WttHk) - 1 (37) 

W4-k - 1 + i) = W'„(«d) - Wy>(«D) + (-)“_^ Vi;;\k) - 1 (38) 


The continuation of double etc. sums is done in a similar manner. In the inner summation one 
uses the continuation of Sa via ■^-function or (^). The results are similar to (^ ) - (|38|) e.g. 

Va,l{-k -1 + S) = Va,l{00) - Kp (OO) 

+ (-)“7rctg(7r(5) ^ITiZiJi(fc) + W^-^\k) - 7rctg(7r5)ITizf (fc)^ - V’(<^)^+ 1e 

In the same manner functions with 2 arguments are continued. There is one shortcoming, 
however, in such a procedure — some of the inhnite sums are divergent, e.g. Wa{oo). This 
may cause problems in some cases. One way out may be that these functions are evaluated as 
hypergeometric series with the corresponding argument, like e.g. 

= r(l + 5)r^2 + 5) ^ 1 + 2 + <5,2 + <5;4). (40) 

Moreover, some of them (e.g. Si) cancel after applying the SW-transformation. 

As an example of this technique we transform the series z^Vs/n into a l/z"" series. Utilizing 
formula (l37|) we get 

^ 1 log^(- 2 ;) + 2 C 3 log(-^) - 3(4 - 2U4 (cx)) + 2 U 3 ,i(oo) - 2 ^ 3,1 (oo) 

n=l ^ 

OO 1 1 / 1 1 -I 

- E - (^Q log'(-^) + ^ logY-z) + {:^ + W 2 {n)) log(- 2 :) 

+ 2(3 + 3IU3(n) + 2114,1 (’^) ~ 2114 ,i(n) H-^ . (41) 

In this formula the hrst line is the residue at n = 0 while the summation is due to residua at 
negative n. 
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8 Conclusion 


For a certain class of selfenergy and vertex functions with only one non-zero mass we have 
developed a new technique of obtaining analytic results. This method uses information from 
the large mass expansion of the diagrams by reproducing the obtained expansion coefficients in 
terms of a certain set of ‘basis elements’, which are harmonic sums in many cases but often more 
complicated (‘W- and V-sums’). However, we did not investigate systematically all 2m-vertex 
functions. In some cases we failed to find the solution. Supposedly some of them include new 
exceptional basis elements. 

Given the expansion coefficients analytically, it is possible to sum most of the series in 
closed form to yield a representation of the diagrams under consideration in terms of poly¬ 
logarithms. Furthermore the analytic form of the expansion coefficients allows by means of 
the Sommerfeld-Watson transformation a direct transition to the large g^-expansion. For di¬ 
agrams involving elliptic functions (essentially those with 3m-cuts) the basis elements of the 
corresponding expansions were not found. 
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Appendix 


A Basic functions 

Here we list all basic functions we used and point their ’transcendentality level’ (TL). 

A.l Zero-fold sums 


1 


(42) 

(43) 


n' 



(44) 


(46) 


where (a = C(®) is the Riemann (^-function. and (^31) have TL= a. (^) and (^) have 


TL= 0. 


A.2 One-fold sums 



( 46 ) 
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( 47 ) 




''■“’■sBJ'? 


All have TL= a. An exceptional sum occurs in the diagram /123 with TL= 2 


Un) = Z 


'2j\ ‘ 1 


^iVi/ jin + l-j) 


A .3 Two-fold sums 


n 1 

j=lJ 

^ (-V 

j=i J 

HVj(n) = - 1). 

= -1) ■ 

If instead of Sa{j — 1) there is an insertion of Sa(2j — 1) we use the tilde, e.g. 


.\-i 


y^An) = t (^/) Y^S,i2j - 1 ), 


etc. All have TL= a + b. 


A .4 Three-fold sums 


Sa,{b+c)in)=Y.—Sb{j - l)Sc{j - 1), 
j=iJ 

Wa,(b+c){n) ^Sb{j - l)Sc{j - 1), 

j=l \ J / J 


All have TL= a + h + c. 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


(56) 

(57) 
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A .5 Functions with two arguments 

Functions S and V may occnr with two argnments 


Sa{n;p) = ^^ 
j = l J 


n 


ya{n]p)=Y^ 

i=i 


2 j\ 

j) j° 


(58) 

(59) 


All have TL— a. 


B The 2-point functions 

Below we present results for the diagrams shown in Fig.2 All fnnctions have argument n — 1 
which is omitted. Fnnctions with two arguments have their hrst argument n — 1 which is 
omitted. Si = Si{2n — 1) and V 2 in /123 is given by (^. 2 ; = /mS. 


1 00 1 

- 
n 


r 


-I 00 

h = -Y.^" 

9 SI 


n 


- log^(— 2 ;)-log(— 2 ;) + C 2 + 2 S '2 — 2 - 1 —- 

2 n n 

2 1 f—V 

- \og^{-z) + - log(-z) - 2 C 2 + 4 iF 2 - 2 - - 2 ^- 


1 00 1 

T .S n‘ 


n \ n 


v-l 


1 00 1 

a3=4e-- 

1 00 1 

/i4 = 4E^ - 
r ^ 


2C2 + 2^2 — 2 


-2 log(-; 2 ) - 3PFi + - 




n 


-Si{;2)\og{-z) 


+2C2 ~ 6V2 + 2V2(; 2) + S'2(; 2) — Ai i(; 2) + 2-1-^^ 

n n 


1 00 1 


r 


n 


1 2 /2n 

-- log(-z) - C 2 + 31/2 + — + 4 
n 


v-l 


n 


1 00 1 

4i23 = ~ E 
r n 


C ,2 ~ 3V2 + 3V2 ) ~ 


f 2 n 


n 




-1 00 

4i25 = ~ E ^"’“2 

g 2 V? 

1 00 1 

4 i 234 = ~ Z^ — 


r 


n 


F2(;4) + 


3 IU 

n 


, n n\ n 


(2n\ ^ 

/ 5 3' 

(J 

- log(-^) + - 

n 


(60) 

(61) 

(62) 

(63) 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


First seven coefficients of the expansion for Ji, Is, J 12 , / 14 , /15 and /125 were compared with 

HI. 
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C Integral representations for the 2-point functions 


Here we give the representations for the 2 -point functions. Following we introduce the 
notation 

c { \ - (-1)“+^ log“(t) log^(l - zt) 


iTi / dt. 

a\ b\ Jo t 

All polylogarithmic functions are particular cases of S functions, namely 

Ua{z) = Sa-l,l{z) . 

We introduce the following new variable {z = /mJ) 


(69) 


(70) 


V = 


l-^^/(^-4) 


1 + JzKz - 4) 


(71) 


Then 


■ Ji = -- log^(-^) log(l - z )-2 log(-z)Li2(2:) 3113(2;) - 6S'i,2(^) 


— log(l — 2; 


( 72 ) 


■ h = 2(2 log(l + z) + 2 log(-2;)Li2(-^) + log^(-;2) log(l + z) + A log(l -b ^)Li2(2;) 

- 2Li3(-z) - 2Li3(z) + 2 Ai,2(z') - 4Ai, 2(;^) - 4Ai,2(-^), (73) 


q^ ■ I12 = Lig {z) - 6C3 - C2 log g - - log^ y-Alogy Lis {y) 

0 

+ 4Li3(?/) - 3Li3(-|/) -b ^Li3(-j/^), 


■ ho = -6S'i, 2(^) - 2 log(l - 2; 


( 74 ) 

( 75 ) 


q^ ■ /i4 = log(2 - 2; 


2;) - 2 log(-2:) log(l - z) - 2112(2: 


- 2 log^(l - h+ log(-^) log^(l - z)- 2(2 log(l - z) 

~ >S'i^2(l/(l ~ ^)^)+2*S'i_2(l/(1 ~ ^))+2*S'i_2( ~ 1/(1 ~ ^)) 
+ ^ log^ y + log^ y (2 log(l + y ^)-3 log(l -y + y^)^ 

- 6C3 - Usi-y^) + ‘^Lisi-y^) - eLk^-y) 


+ 2 log2/ fLi2(-|/^) - U2{-y^) + 3Li2(- 


( 76 ) 


q^ ■ Ii5 = 2Li3(^) - log(-2;) 112(2;) -b C2 log(l - z) 


1, 3 1, 2 

g log y - - log y 


y) - 3 \og{l - y + y^)J 
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( 77 ) 


- 6C3 - ^Li 3 (-|/^) + 3Li3(-|/) + 8 Li 3 (|/) 

+ logy (^Li2(-|/^) - 3 U 2 {-y) - 8Li2(?/)^ , 

■ h23 = -C 2 (log(l -z)+ log y'^ - 6(3 - ^ log(l - y + y^) log^ y 

+ ^h{-y^) - 9Li3(-?/) - 2 log?/ (^Li2(-?/^) - 3Li2(-?/)^ , ( 78 ) 

■ /i 25 =-2 log^ ?/ log(l -y)- 6C3 + 6Li3(?/) - 61 og?/Li 2 (?/), ( 79 ) 

q^ ■ ^1234 = -6C3 - 12Li3(?/) - 2 AU^{-y) 

+ 8 log?/ (^Li2(?/) + 2Li2(-?/)^ + 2 log^ y (^log(l - ?/) + 2 log(l + y)^ . ( 80 ) 

Different analytic representations for Ji, J13, /125 and /1234 were obtained in p[. In [Q 
the same four integrals were calculated for arbitrary dimension d in terms of hypergeometric 
functions. and J12 were calculated also in [^]. The numerical check shows agreement for 
these integrals. 


D 3-point functions 

Below we present results for the diagrams shown in Fig .3 All functions have argument n — 1 
which is omitted. Overlined functions have argument 2 n — 1 which is omitted (S'! = S'i( 2 rr — 1 ) 
etc.). 14,1 in P56 is given by (|^) and z = 


Pi = 


(q 


2\2 


E. 

77=1 


n 


-+ -S 2 - -Si H— -log(- 2 ;) + - log^( 

£ L 2 2 2 n 2 


+ (C 2 - 4^2 - 2 E _ log(-z) + (Si + log^(- 2 ;) - ^ log^(-z)| 


-1 00 

F3 = -4.E-" 


(q 


2\2 


n=l 


Irl 
+ £L 2 


3 


( 81 ) 


+C4 + C 3 F 1 + ^( 2^2 + C2S't + + 3Fi,3 - 2Fi,(1+2) - 8 F 1 F 3 - 2FiFi,2 - 2 F 2 " + 2S'tS2 


+ 




(—C3 ~ 2C2S'i + ^3 + 2 S'i 52 ) log(—^) + ^(^2 — 2*^2) log^(—2)| 

_)n 


( 82 ) 


(q 


2\2 


E. 

77=1 


n 
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i ~6(^3 + 2(^2>S'i + QS^ — 2S'iS'2 + 4-^ + 2— 

[ n n 

+ (~4S'2 + — 2—) log(—+ Si log^(—^)| 


( 83 ) 


Pfi = 


E 


(g 2)2 n 

{4h- + iog(“^) 

(2 — 3S'2 + 4-^2 ~ 3- - + f3S'i H—^ log(— 2 :) — - log^(— 2 :) 

/Z TX Tlj ^ 

C 7 9 

+ 2((^3 + 1 C, 2 S\ + 2 -25*3 — 9 S'iS '2 + + 12 ii' 2 ,i + 4 S'i-K '2 ~ 7 : — 

n in 

QS! Si 7 , / _ , 7 ^ , 9^2 , .^1 7 




~2S ■ ■ 75 + + 2 ^^- + 2 ^"+ 


Pl3 = 


1 


(iY s 

f 11 


s—7-P2 + 4--P3-2P12 

I ^2 2 ^ e L 2 

+42*^2 - ^^4 - Pi, 3 + 3Pi,(i+2) - P1P3 - 8PiPi,2 + ^^2' + ^P?P2 


Pl2 = 


E 


'2n\ ^ 1 




2 1 r 

£2 e L 


-614^1 + 2Pi H-2 log(— 2 ;) — 614^2 — 1844^1 ,1 — 13^2 


n 


+Pi — 6Pi44^i + 2 - 1 —- -\- 2 Pi-^ log(—+ log^(—^)| 

Tt ‘Th T~t I 


T 00 /_ \n 

^6 = — 


Si " 

—6((^2pi + 64/3 — 2P3 — 10773 — 124/2^1 + 124^ 1 — 12772^1 + 12Pi772 + 6 


4P 


n 


S 2 


3 ^ - 4 E + (^-QV 2 + 3 P 2 + 4^) log(-s) - 2Pi log2(-s) 


E' 

n 


Pl 26 — 


'2n\ ^ 1 


00 


—^ H— 

^2 £ L 


15 


-Pi - log(-s) I - ^P 2 - y p? + 4PiPi 


+2— - Pi log(-s) + ^ log^(-^) 
n 2 


(84) 


(85) 


( 86 ) 


(87) 


( 88 ) 
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75*3 — 4 >S’i ^2 — ST^a + 8S1K2 + 85*3 + 86*1^2 + 87^2,1 + 4S’iS’2 — 8S1K2 — 8S1S2 


+ (-2^2 + 47^2) log(-;^)| 

1 00 

(9 ) 

I ^ s, - i log(-2) + i -ife + 3S, + 2Sl - 2S, log(-2) + i log"(-j) 
+- -ICs - C 2 S 1 + ^Ss + 2S,,2 + lOSiSj + |s? 

o o o 

+ (^C 2 - 4^2 - 45 i) log(- 2 :) + 2^1 log^(-^) - ^ log^(-^) 

^7 32 ^2 

“16^" - y Ca^i + C2^2 - 2 C 2 SI + 4^4 - 45i,3 - 2^1,(1+2) + 75^ + -S,S^ 

+8S,S,,2 + 17 S^,S 2 + 1^1" + (^Cs + 2(2^1 - ‘^83 - I6S1S2 - log(-;^) 




1 ) log^(-^) - ^-^1 log^(-^) + ^ log^(-^) 


1 00 

^ 3 =T^E^"(-r 

e ~ 2C2 >S'i + 2S3 — 2K3 + 4 S'i 772 + (C2 ~ >S'2) log(— 

— C4 ~ 2(381 — 7(282 — 4C2 *S'i + 7^27^2 
—84 + ^-8i + Q8183 + 281,3 + 8774 - 881K3 + 88fK2 

+ (C 3 + 4 C 2 >S'i — 81,2 — 88382 — 4773^ log( —+ (^—(2 + -82 + K2) log^( —^)| 

1 00 1 

y e L 2 2 n J 

- 5 C 2 ^i + ^^3 + 2 ^ 1,2 - ^ 5 i ^2 - ^^ 8 l + :^ + 7 ^ - 4^1 

6 2 b n n j 

I -4w"i + - I-4W2 - I2W3 1 - 12^114^1 - I652' 

[ e ’ J 

— 18(2Wi — 4IV3 — 85*3 + 1214^12 — 844^2,1 — 1244^1 (i+i) 
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- 16 ^ 1,2 - 8 ^ 11^2 - 24 ^ 11 ^ 1,1 - 16 ^ 1^2 - 1251 ^ 1^1 


( 93 ) 


E Useful sums 


In Appendix B we have presented our results for two-point functions in terms of the basis 
elements given in appendix A. Summing the series, we were able to obtain closed expressions 
in terms of polylogarithms. In appendix D we also gave the series representation for the vertex 
functions under consideration. The summation and representation in terms of integrals of 
polylogarithms is done for separate structures, occuring repeatedly in these series. 

In the case when only harmonic sums enter the expression the formulae given in |^] are 
very useful. For convenience we introduce the notation 


Lii(z) = -log(l - z). 


( 94 ) 


E.l Sums involving S and K 

The following general relations hold true 




n=l 

oo 




10 Z' —, n 


i:/(n - 1 )—i: 


n=l 


n=l 


f{n) -/(n- 1) 




where f{n — 1) is any function. 
As a particular result we have 


OO 1 


n=l 

oo 




Sain - 1) z*" = Lia(z), 

n=l ^ 


n=l 

oo 




Kain -l)z^ = Lia(- 2 ;) 

n=l ^ 


(95) 

(96) 


(97) 

(98) 

(99) 
( 100 ) 


Only S'-functions (for brevity: on the l.h.s. Sa = S'a(n—1) etc. and on the r.h.s. LF = lAaiz) 
etc.) 


OO c 

^^^- = LiiLi2-2^1,2, 

OO o 2 1 

^ — z” = LiiLi2 -|- xLii — 2S'i 2 , 

OO q 1 

.3 " 2 


( 101 ) 

( 102 ) 

(103) 
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/ , -^ — ■;^Li]^Li2 + L11L13 — -L12 — Lii^i 2, 

TX Zi Z 


n=l 
00 ^3 


^ = llit + LiiLis - + hijU^ - 3Lii5i,2, 

n=i n 4: z z 


“ ^1.2 . 1 . 


E — = oLi?Li2 - 2Lii5i,2 + 3 ^ 1 , 3 , 

00 Q -I 

E|-” = ^Li^- 2 S., 2 . 

n=\ ^ 

00 c 2 1 

i:; 5 "” = 5 Li? + 2 Sw- 2 S: 


n=l 

00 


n 


E -s? *-" = 




n=l 

OQ 


Y.SiS2 Z^ = 

n=l 

00 

EE 2 ^" = 

n=l 

00 

Y.s\z’' = 


l-z' 

z 


l-z 

-^1 

1 -; 2 \ 

/ 
z 


2,2, 

Lii + Li 2 ^, 

^Lis + LiiLi2 — 5'i^2^ , 
LiiLi 2 — 2 S'i 2 ^ , 


n=l 

00 


J2s,s,z^ = 


l-z' 
z 


Li3 + 3LiiLi2 + Li^ — 35*12^, 
1 . 


n=l 

00 


l-z 


Li 4 + LiiLi 3 — *^ 2,2 


Y.sh"^ 


n=l 

00 


l-z 


Li 2 + Li4 — 45*2,2 1 5 


5*^5*2 j [ Li^Li2 + Li4 + 2LiiLi3 — 2LiiS*i^2 + 2S*i^3 ~ 25*2,2) ? 


n=l 

00 


(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 
( 112 ) 

(113) 

(114) 

(115) 


Y.Stz" = 


z 


n=l 

00 


ESi,3^” = 

n=l 

00 

n=l 

S and K 
^2 


l-z 

z 


Li2 + Li4 + 4LiiLi3 + 6Li^Li2 + Li^ — 12Lii5*i^2 + 125*1^3 ~ 85*2, 2 J ^(IIS) 
1 , 


l-z 


LiiLi 3 — -Li 2 1 , 


= 


^ ' - • - • ; 2 t ; ;2 


l-z 


LiiLis + -Li3Li2 — ~ LiiS' 1^2 


“ K- 




n=l 

00 


n 




n=l 

00 


n 


E —= 

^ n 


n=l 

“ SiK. 


Jo 1 — t 

Li 3 (-t) 
/o 1 — t 


dt, 


dt, 


E z^ = - log(l - z)(j){z) - [ 


log(l — t)Li2(—t) + Li3(—t) + >S'i,2(~^) 
1-t 


(117) 

(118) 

(119) 

( 120 ) 
( 121 ) 

(ft, (122) 
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^ K. 


= 0 ( 2 :) log ^ + 


logtLi2(-t) 
'0 1 — t 


f; S,K2 z- = - S,,2{-z) - Li3(-;^)'), 

n=l ^ V / 

f; S,K, z- = ^ (-U,{-z) - S2,2{-z) - r dt], 

1 - ^ V Jo I-t J 

°° z ( 1 

^ S2K2 z^ = -U,{-z) - -Lili-z) + 2 S2,2 {-z) + <p{z) \ogz 

n=l ^ V ^ 


logtLi 2 (-t) 
/o 1 — t 


-^ 3,1 z"^ - ^ ^ |^-*S'2,2(-2:) j , 


(123) 


(124) 


(125) 


(126) 


(127) 


(128) 


00 z ( 1 

SIK 2 z"^ = -— -Lu{-z) - -Lili-z) - 2 Si^ 3 {-z) + 0 (^) (log z + 2 log(l - z)) 

n=l ^ \ ^ 


where 


(2 log(l -t) - logt)Li 2 (-t) + 2 Li 3 (-t) + 2 S'i, 3 (-t) 


4>{z) = -Si^ 2 {z‘^) - Si^ 2 {z) - Si^ 2 {-z) + log(l - ^)Li 2 (- 2 :). 


(129) 


(130) 


E.2 Sums involving V and W 

The factor can be written in the following useful form 


^2n\ n d 

, n ] 2 Jo 


dss^-\l-sy 

Jo 


Again we make use of some general formulae. 


For any function / we have (hereafter p = s(l — s)) 


2 n\ /(n - 1 ) 


— 1 'j /*! r/s 

/ -Yfin-I){zp)\ 
n Jo s 


£ (£1 h' - *■”=T ("p)” ■ 

n=l \ 'V Z JO b 

“ /2nYV(2n-l) „ 1 A ds “ ., . r , 1 

i:(,j =-J„ vi:/("-i)[(v^) +(-vzp)] 

When the function Va^iji — 1) contribute we have 


^^2 n [^ds 

Y—^= —h{z,p) 

“A n Jo s 


( 131 ) 
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°° 1 r - 1 /•! r/s 

V 2 + V 2 z'^ = j — log(l - z) log(l - pz) 

J Jo s 

- I2{z,p)\0g{l-p)+l2,l{z,p) 

n Jo s V J 

^ - V" 2 ,l - 14,1 = 7 ; — log(p)/ 2 (^,p) 

^inL J 2Jo s 

^ V ds \ 

XI — 2 :” = - / — log(l - z)l.i2{zp) + 5'i,2(2:) + Si^2{zp) + 13 ( 2 :,p) 
n Jo s L J 

00 /*! /7 q r 

/ — ^2(^,p)log(2:)+log(l-p)(Li2(l-^)-C2) 

n=l ^ JO S L 

- ^log^(l -P)l0g(l -^) +4l(^>P) +^1,2(^,P) , 

where 


h{z,p) = U 2 {^——) - U 2 {p) - log(^— —)\og{l -p) + ^log 2 (l - zp) 
1 — zp 1 — z 2 

h{z,p) = Li3(p) - C 3 - Li3(4 ^——) + Li 3 ( 7 —-) 

1 — zp 1 — zp 


+ log(^- -) [Li 2 (^j- —) - Li 2 (^-^)] - - log(^) log 2 (-- -) 

1 — zp i 1 — zp 1 — zp \ 2 1 — zp 

fr){ 1 _ 2''\ 1 1 _ 2J, 1 
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1 — zp 2 1 — zp 2 
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- - log (1 - zp) 

0 
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1 — zp 

111 
+ o log(l - zp) log^(^(l -p))-77 log(^(l - p)) log^(l - + 7 log^(l - ^p) 

2 2b 


When the functions Wa,b{n — 1) contribute we have 


3 X “^1 z'^ = - j ds—^ — log(l - zs) 

\ n J n Jo 1 — zp 

3XP^1 -SiWiz"^ = -[ dsz^ —log(l- 2 ;s)log('— 
J n Jo I - zp V (1 

°° f2v\~^ 1 r 1 /•! z r 

X ( ) ~ 3hhi 1 + hh 2 z^ = — ds- - 2 Li 2 (p 2 :) + 

^ \ n J ni J Jo 1 — zpi 
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n n 


1 /*! ^7 9 r 9 

—Wi z” = - / — Li 2 (s) - Li 2 (--) + log(l - s) log(l - zp) 

Jo p l 1 — zp 


^log^(l - z. 
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log(l — zp) log(l — s) log 


°° '2n\ ^ 1 


E 

n=l 



n 


31^11 + W 2 


z = 


1-sJ] 

1 '^ ds 
0 p 


“ ^2n\ ^ 1 


“ 2n\ h 


- log 

E 

n=l 
c 

3E 
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zs 


1 — zp \{1 — s){l — zp) 


+ 2U3{pz) + 5 ' i , 2 ( p 2 :) 


5 

+ 7 log^(l - ^P) 

0 


9$5(^,p) = 2 

-4 


Li 3 (p 2 :) - 65'i,2(p 2:) - 3 log(l - 2 :p)Li 2 (p 2 :) - - log^(l - zp) 

0 


zs‘^ zs‘^ 

Li3(T^)+2^l,2(: 


1 — zp 


Li 


zs 


1 - zp^\ - zp' ^^V(l - s)(l - zp)^J • 


24 



References 

[1] G. Kallen and A. Sabry, Dan.Mat.Fys.Medd. 29 (1955) No. 17; R. Barbieri, J.A. Mignaco 
and E. Remiddi, Nuov.Cim. All (1972) 824, 865. 

[2] D.J. Broadhurst, Z.Phys. C47 (1990) 115. 

[3] D.J. Broadhurst, J. Fleischer and O.V. Tarasov, Z.Phys. C60 (1993) 287. 

[4] S. Bauberger et al., Nucl.Phys. B434 (1995) 383; V. Borodulin and G. Jikia, Phys.Lett. 
B391 (1997) 434. 

[5] R. Scharf and J.B. Tausk, Nucl.Phys. B412 (1994) 523. 

[6] A. Ghinculov and J.J. van der Bij, Nucl.Phys. B436 (1995) 30; V. Borodulin, G. Jikia, 
Phys.Lett. B391 (1997) 434. 

[7] A.V. Kotikov, Phys. Lett. B254 (1991) 158; B259 (1991) 314; B267 (1991) 123 (Err. 
B295 (1992) 409); In: Artificial Intelligence and Expert Systems for High Energy and 
Nuclear Physics, (Oberamergau, Germany, 1993) p.453; JINR preprint E2-93-414. 

[8] J. Fleischer, A.V. Kotikov and O.L. Veretin, Phys.Lett. B417 (1998) 163. 

[9] D.J. Broadhurst, |hep-th/9803091 . 

[10] F.V. Tkachov, Preprint INR P-0332, Moscow (1983); P-0358, Moscow 1984; 

K.G. Ghetyrkin, Teor. Math. Phys. 75 (1988) 26; ibid 76 (1988) 207; Preprint, MPI- 
PAE/PTh-13/91, Munich (1991); V.A. Smirnov, Gomm. Math. Phys. 134 (1990) 109; 
Renormalization and Asymptotic Expansions, Birkhauser, Basel, 1991. 

[11] J. Fleischer and O.V. Tarasov, Z.Phys., C 64 (1994) 413; J. Fleischer, V. A. Smirnov and 
O. V. Tarasov, Z.Phys. C74 (1997) 379; J.Fleischer et al, Eur.Phys.J. G2 (1998) 747. 

[12] A.I. Davydychev and J.B. Tausk, Nucl.Phys. B397 (1993) 123. 

[13] A.I. Davydychev, V.A. Smirnov and J.B. Tausk, Nucl. Phys.B410 (1993) 325; F.A. 
Berends, A.I. Davydychev and V.A. Smirnov, Nucl. Phys.B478 (1996) 59; L.V. Avdeev 
and M.Yu. Kalmykov, Nucl. Phys.B502 (1997) 419; G. Weiglein, |hep-ph/9711254| . 

[14] F.A. Berends, A.I. Davydychev, V.A. Smirnov and J.B. Tausk, Nucl. Phys.B439 (1995) 
536; 

[15] F.J. Yndurain, Quantum Chromodynamics (Springer-Verlag New York Inc., 1983). 

[16] L. Lewin, Polylogarithm and associated functions (North-Holland, Amsterdam, 1981); 

N. Nielsen, Nova Acta Leopold. 90 (1909) 123. 

[17] J.A.M. Vermaseren, NIKHEF-98-14, Piep-ph/9806280| . 

[18] S. Bauberger and M. Bohm, Nucl.Phys. B445 (1995) 25. 

[19] D.I. Kazakov and A.V. Kotikov, Nucl.Phys. B307 (1988) 721; Theor.Math.Phys. 73 (1988) 
1264; A.V. Kotikov, Theor.Math.Phys. 78 (1989) 134. 


25 









[20] D.I. Kazakov and A.V. Kotikov, Phys.Lett. B291 (1992) 171. 

[21] E. Remiddi, Nuovo Cim. IlOA (1997) 1435. 

[22] A.V. Kotikov, |hep-ph/9807440| . 

[23] J.A.M. Vermaseren, Symbolic Manipulation with FORM (CAN, Amsterdam, The Nether- 
lads, 1991). 

[24] A. Devoto and D. W. Duke, La Rivista del Nuovo Cimento, Vol. 7, No. 6, pg.l (1984). 

[25] K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159; F.V. Tkachov, Phys. 
Lett. BlOO (1981) 65; A.N. Vassiliev, Y.M. Pis’mak and Y.P. Honkonen, Theor. Math. 
Phys. 47 (1981) 465. 

[26] V.A. Smirnov, Phys. Lett. B404 (1997) 101; 


26 





Figure captions 

Fig. 1 : 2-loop selfenergy- and vertex diagrams considered in the paper. Line numbering and 
kinematics. 


Fig. 2 : 2-loop selfenergy diagrams evaluated in this work. Solid lines denote propagators 
with the mass M; dashed lines denote massless propagators. 


Fig. 3 ; 2-loop vertex diagrams evaluated in this work. Kinematics as in Fig. 1 = P 2 = 0). 

Solid lines denote propagators with mass M; dashed lines denote massless propagators. 







Fig. 3 





Fig. 2 



